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ABSTRACT 

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter 
HD 189733b obtained in the 3.6 and 4.5 /iin bands using the Spitzer Space Telescope. When combined 
■ with previous phase curve observations at 8.0 and 24 /im, these data allow us to characterize the 

exoplanet's emission spectrum as a function of planetary longitude and to search for local variations in 
its vertical thermal profile and atmospheric composition. We utilize an improved method for removing 
the effects of intrapixel sensitivity variations and robustly extracting phase curve signals from these 
data, and we calculate our best-fit parameters and uncertainties using a wavelet-based Markov Chain 
Monte Carlo analysis that accounts for the presence of time-correlated noise in our data. We measure 
a phase curve amplitude of 0.1242% ± 0.0061% in the 3.6 (im band and 0.0982% ± 0.0089% in the 
4.5 /im band, corresponding to brightness temperature contrasts of 503 ± 21 K and 264 ± 24 K, 
P I 1 respectively. We find that the times of minimum and maximum flux occur several hours earlier than 

r~r1 . predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas 

by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from 
our previous observations at 8 /im, and we present new evidence indicating that the flux minimum 
observed in the 8 /im is likely caused by an over-shooting effect in the 8 /im array. We obtain improved 
q ■ estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% ± 0.0040% in the 3.6 /im band 

?-h " and 0.1787% ± 0.0038% in the 4.5 /im band, corresponding to brightness temperatures of 1328 ± 11 K 

\ and 1192 ±9 K, respectively; these are the most accurate secondary eclipse depths obtained to date 

for an extrasolar planet. We compare our new days ide and nightside spe ctra for HD 189733b to 
the predictions of ID radiative transfer models from iBurrows et al.1 (|2008[ L and conclude that fits 
to this planet's dayside spectrum provide a reasonably accurate estimate of the amount of energy 
transported to the night side. Our 3.6 and 4.5 /im phase curves are generally in good ag reement with 
the predictions of general circulation models for this planet from lShowman et al.l (|2009|) . although we 
require either excess drag or slower rotation rates in order to match the locations of the measured 
qq . maxima and minima in the 4.5, 8.0, and 24 /im bands. We find that HD 189733b's 4.5 /im nightside 

flux is 3.3er smaller than predicted by these models, which assume that the chemistry is in local 
thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which 
should lead to an excess of CO and correspondingly enhanced 4.5 /im absorption in this region. This 
result is consistent with our constraints on the planet's transmission spectrum, which also suggest 
excess absorption in the 4.5 /im band at the day-night terminator. 
Subject headings: binaries: eclipsing — planetary systems — techniques: photometric 
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1. INTRODUCTION 

Observations of eclipsing extrasolar planetary systems, 
in which the planet periodically passes in front of and 
then behind its host star, have proven to be a powerful 
diagnostic tool for studies of exoplanetary atmospheres. 
Because the probability of a transiting geometry scales 
as R±/a where i?* is the radius of the host star and a 
is the planet's semi-major axis, the majority of currently 
known transiting planet systems have orbital periods of 
just a few days. At these distances, the time scale for the 
planet to achieve synchronous rotation is short compared 
to the typical ages of the systems, leading to the predic- 
tion that a maj ority of these transiting p lanets should be 
tidally locked (|Bodenheimer et al.1 120011 ) . In this paper 
we focus on the class of gas giant planets known as "hot 
Jupitcrs" , which typically have orbital periods on the or- 
der of 1 — 3 days and atmospheric temperatures ranging 
between 1000 - 3000 K. 

One fundamental question for these planets is what 
fraction of the incident flux absorbed on the planet's 
day side is subsequently transported around to the 
night side. Atmospheric circulation models predict 
that these planets should develop a broad superrotat- 
ing (eastward) equatorial jet that circulates energy be- 
tween the day and night sid e s (e.g., [S howman & Guillot 
2002t iShowman et all 1 120091 ; IShowman fc Polvanil 1 20111: 
Langton fe Lauehlinl [20081: iPobbs-Pixon et all I201C 
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2009). Several recent ground-based studies 
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2010; Waldm amTet al.ll2012D have also re- 
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Heng et all 120111: iRauscher fc Menoul l2012b| ). Pepend- 



ing on the relative strengths of these winds, these plan- 
ets could exhibit large gradients in both temperature and 
composition between the two hemispheres. We can con- 
strain the efficiency of the day-night circulation by mea- 
suring changes in the infrared brightness of the planet 
as a function of orbital phase; the day-night brightness 
contrast can then be translated into a day-night temper- 
ature contrast. There are currently well-characterized 
phase cu rves published for seven planets, includin g 
v And b (Harrington et all 120061: iCrossfield et al]l25l0h . 
HP 189733b (IKnutson et al.l 120071 l20 09a|V HP 149026b 
Knutson et al.l l2009cJl . HP 806 6b iLauehlin et al " 



HAT-P-7 b~(IBorucki et all I2009L IWelsh et al 



20ll: CoRoT-lb ([Snellen et al.1 12009h~~and WASP- 12b 
(|Cowan et aLll2012[ ). with" more sparsely sampled phase 
curves for three additional plan ets (51 Peg b, HP 
209458b, and HP 179949b) from ICowan et all (pOOl . 
These data indicate that hot Jupiters display a diver- 
sity of circulation patterns, ranging from relatively small 
day- night temperature gradients (e.g., HP 189733b) to 
large temperature gradients (e.g., WASP- 12b). 

Of these four systems, HP 189733b stands out both 
as having the best-characterized phase variation, and 
also as the only system with phase curve observations 
at more than one wavelength. We know more about 
this planet's atmosphere than that of any other extra- 
solar planet; ke y results include t he detection of a high- 
altitude haze dPont et al.l 120081: ISing et all 1201 lft and 
sodiu m absorption (jRedfield et al.1 120081 : iHuitson et al.1 
2012) in its visible-light transmission spectrum, as well 
as more controversial detections of methane, carbon 
monoxide, and wa ter absorption in i t s infrared trans- 
mission spectrum dSwain et al.l l200"8t ISing et al.l 120091 : 
iPesert et al.l 120091 : iGibson et al.l l2011al lbl) and carbon 
dioxide absorption in its dayside emission spectrum 



ported detections of methane emission from the planet's 
day side, alth ough these results ha ve been the subject of 
some debate (|Mandell et al.l l20l"lT ) . HP 189733b's day- 
side emission spectrum has also been characterized in the 
near- and mid-infrared using bot h Spitzer IRAC photom- 
etry dCharbonneau et al.l |2008|) and IRS spectroscopy 
(|Grillmair et al.l 12008). with additional constraints on 
its var iability in the 8 /im IRAC band from lAgol et all 
(|2?)T0l) . 

Pespite the extent of the data available for this planet, 
there are still a number of open questions regarding the 
properties of its atmosphere. The single largest out- 
standing question centers on the issu e of whether or not 
the chemistry is in equ ilibrium (e.g.. iMoses et al.l 120111 : 
iVisscher"^ ; Moses 2011); it is likely that the atmospheric 
circulatio n plays an important role in shaping this chem- 
istry (e.g.. lCooper fc Showm an 2006). Although we have 
observational constraints on relative abundances for the 
planet's day side and the day-night terminator, we know 
very little about the properties of HP 189733b's night 
side. In this paper we present new full-orbit phase curve 
observations for HP 189733b in the 3.6 and 4.5 /im bands 
obtained with the Spitzer Space Telescope during its ex- 
tended warm mission. We combine these data with pre- 
vious observations at 8.0 and 24 jim. to provide the first 
detailed characterization of its emission spectrum as a 
function of orbital phase. Our data include two sec- 
ondary eclipses and one transit in each band, which we 
use to derive improved estimates for the planet's orbital 
ephemcris and dayside emission spectrum. 

2. OBSERVATIONS 

We analyze two full-orbit phase curves for HP 189733b, 
obtained in the 3.6 and 4.5 mi cron bands using the IRAC 
instru ment dFazio et al.ll2004f) on the Spitzer Space Tele- 
scope (W erner et al.ll2004D . The 3.6 /im observation was 
obtained on UT 2010 Pec 27 - 29 and has a total du- 
ration of 68.5 hours, corresponding to 1,722,624 images. 
The 4.5 /Ltm observation was obtained on UT 2009 Pec 
22 — 24 and has a total duration of 66.7 hours, corre- 
sponding to 1,711,872 images. Pata in both bandpasses 
were obtained in subarray mode with 0.1 s exposures in 
order to avoid saturation and to minimize the total data 
volume. As a result of the limited on-board memory, 
both observations required three breaks for downlinks 
with a duration of 1 — 2 hours per downlink. In order 
to downlink data the spacecraft must point towards the 
earth and then reacquire the target, which produced off- 
sets of up to half a pixel in the star's position after each 
downlink (Fig. [T]and[2|). 

We calculate the B JP_UTC values at mid-exposure for 
each image using the MBJP_OBS keyword in the image 
headers. Each set of 64 images obtained in subarray 
mode comes as a single FITS file with a time stamp cor- 
responding to the start of the first image; we calculate 
the time stamps for individual images assuming uniform 
spacing and using the difference between the AINTBEG 
and ATIMEENP head ers, which record the s tart and end 
of the 64- image series. lEastman et~aT1 (|2010l ) further ad- 
vocate a conversion from UTC to TT timing standards, 
which includes a more consistent treatment of leap sec- 
onds. We note that for the dates spanned by these obser- 
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Figure 1. Measured x positions (top panel), y positions (upper 
middle panel), and noise pixel values (lower middle panel) as a 
function of orbital phase for the 3.6 fim phase curve observations. 
Data has been binned into two minute intervals. The raw pho- 
tometry is shown in the bottom panel. Gaps in data are due to 
spacecraft downlinks, and the offset in positions between downlinks 
is a result of the need to re-acquire the star after each downlink. 

vations the conversion from BJDJJTC to BJD_TT sim- 
ply requires the addition of 66.184 s, and we proceed us- 
ing BJDJJTC in order to ensure consistency with other 
studies of this planet. 

2.1. Photometry 

Subarray images have dimensions of 32 x 32 pixels, 
making it challenging to estimate the sky background 
independent of contamination from the wings of the 
star's point spread function. We choose to exclude pix- 
els within a radius of 12 pixels of the star's position, as 
well as the 13th-16th rows and the 14th-15th columns, 
where the stellar point spread function extends close to 
the edge of the array. We also exclude the top (32nd) row 
of pixels, which have values that are consistently lower 
than those for the rest of the array. Lastly, we mask 
out a 4 x 3 pixel box centere d on the position o f HD 
189733b's M star companion (|Bakos et al.ll2006aD . We 
then iteratively trim 3er outliers from the remaining sub- 
set of pixels, create a histogram of the remaining values, 
and fit a Gaussian to this histogram to determine the sky 
background for each image. We find that the background 
contributes 0.2% of the total flux at 3.6 /mi and 0.1% of 
the total flux at 4.5 /mi. 

We determine the position of the star on the ar- 
ray in each image us i ng flux-weighted centroidin g (e.g., 
IKnutson et all 120081: ICharbonneau et al.l I2008D . We 
also tried fits using Gaussian position estimates (e.g., 
iStevenson et al.l 120101: lAgol et al.ll2010f ). but found that 
this method results in an inferior correction for the in- 
trapixel sensitivity variations in our data. Before esti- 
mating the position of the star on the array we first cor- 
rect for transient hot pixels in each set of 64 images by 



Figure 2. Measured x positions (top panel), y positions (upper 
middle panel), and noise pixel values (lower middle panel) as a 
function of orbital phase for the 4.5 /im phase curve observations. 
The raw photometry is shown in the bottom panel; see Fig. [l]for 
a complete description. 

replacing outliers more than 3<r away from the median 
flux at a given pixel position with this median value. We 
then subtract the best-fit background flux from each im- 
age and iteratively calculate the flux-weighted centroid 
for a circular region with a radius of either 4.0 (3.6 /im) 
or 3.5 (4.5 /im) pixels centered on the estimated position 
of the star. Increasing or decreasing the size of this region 
does not significantly alter the time series but does result 
in a slightly higher scatter in the normalized light curve. 
We also estimate the width of the stellar point spread 
function in each image by cal culating a quantity known 
as th e noise pixel parameter ([Mighell 20Q5| ILewis et al.l 
HoH), which is defined in Section 2.2.2 ("IRAC Image 
Quality") of the IRAC instrument handbook as: 

= %tT (1) 

where is the measured intensity in the i th pixel. This 
noise pixel parameter is equal to o ne over the sharpness 
param eter S\ , first introduced by iMuller fc Buffingtonl 
i|1974D . and proportional to th e FWHM of the stellar 
point spread function squared (Mighell 2005). We find 
that we obtain the best results in both bands when we 
calculate f3 using a circular aperture with a radius of 4.0 
pixels. 

We calculate the flux in each image using aperture 
photometry with either a fixed radius ranging between 
2.0-5.0 pixels in 0.1 pixel steps or a time-varying radius 
equal to the square root of the noise pixel parameter 
f3 with either a constant scaling factor ranging between 
0.8 — 1.2 or a constant offset between —0.4 and +0.4 
pixels. We optimize our choice of aperture individually 
for each of the four data segments defined by the down- 
link breaks, as the star falls on a different region of the 
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pixel in each section. We find that for the 3.6 /jm array 
we obtain the best results for apertures with radii equal 
to ^/(/3)+[0.1, 0.2, 0.2, 0.2] pixels respectively in our four 
segments; this corresponds to a range of 2.3 — 2.6 pixels 
with a median value of 2.5 pixels. In the 4.5 /im band we 
prefer a fixed aperture with a radius of 2.3 pixels in all 
segments. We obtain consistent results for the best-fit 
phase curve and eclipse parameters in both bands over a 
range of either fixed (4.5 //m) or time-varying (3.6 /im) 
aperture sizes, albeit with a larger standard deviation 
and correspondingly higher uncertainties for the less op- 
timal apertures. We remove outliers from our final light 
curves using a moving median filter in flux with a width 
of 50 points, where we discard points that lie more than 
5cr away from the median flux values. This corresponds 
to 0.008% of our data at 3.6 /mi and 0.025% at 4.5 mi- 
cron. 

In addition to the intrapixel sensitivity variations de- 
scribed in §2.2[ our data also exhibit a short duration, 
ramp-like behavio r similar to that observed in the IRAC 
8.0 u rn array fe.g. lKnutson et al.ll2007l l2009d: lAgol et al.l 
2010) at the start of each phase curve observation and 
again after each downlink break. We find that the am- 
plitude of this ramp decreases in each successive seg- 
ment, with the largest ramp at the start of the first 
segment. This effect may be due to charge-trapping in 
the a rray, as has been sug gested for longer wavelengths 
(e.g. JKnutson et aLll2007T) . or it may be related to a set- 
tling of the telescope at a new pointing. This effect has 
an asymptotic shape and generally converges to a con- 
stant value on time scales of an hour or less; we elect 
to trim the first 60 minutes of data at the start of the 
3.6 /jm phase curve observation and the first 30 minutes 
of data at the start of the 4.5 /im phase curve observa- 
tion, repeating this trim after each downlink break. As 
before, we select our trim intervals in order to minimize 
the standard deviation of the residuals from the best-fit 
solution. Such trimming is standard practice for Spitzer 
secondary eclipse observations in these bands, which 
typically have a duration of less than eight hours (e.g. 
Knutson et"aHl2009bt iTodoroy et "al~ll2010t iFressin et al.l 
20101 : lO'Donovan et al.l 120101: iDeming et al.l l2011h . We 



also tried including exponential functions at the start of 
each data segment in our fits to the 3.6 phase curve, 
which displayed a stronger ramp than our 4.5 /im data, 
but we found that these functions did not significantly 
improve the quality of the fits, nor did they change the 
best-fit phase curve amplitude and the times of minimum 
and maximum flux. 

2.2. Correction for Intrapixel Sensitivity Variations 

Fluxes measured at these two wavelengths show a 
strong correlation with the changing position of the star 
on the array. This effect is due to a well-documented 
intra-pixel sensitivity var i ation (e.g. 



Reach et al.l 200E 
20051 120081: iMorales-Calderon et al.l 
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2006; Knutson et al 



2008j) , in which the sensitivity of an 



individual pixel varies between the center and the edge. 
We find that quadratic functions of x and y position, 
which are commonly used to correct shorter (< 10 hour) 
transit and secondary eclipse observations, are inade- 
quate to describe the pixel response over the range of 
positions spanned by our phase curve observations. We 



instead correct for these sensitivity variations by approx- 
imating the star as a point source on the array and using 
the measured fluxes as a function of x and y position 
to create a map of the pixel respo nse. Our approach i s 
based on the method described in iBallard et al.l (|2010f ). 
but includes a number of changes that improve the per- 
formance and reduce computational overheads for our 
long light curves (75 hours, vs. approximately 20 hours 
in Ballard et al.). We calculate the effective pixel sensi- 
tivity at a given position as follows: 



xe 



y z - Vj f/2a 2 yJ 



i=0 



x e 



r $7i 

(2) 




where F meas j is the measured flux in the jth image, Fq j 
is the intrinsic flux, Xj, yj, and f3j are the measured x 
position, y position, and noise pixel values, and o~ x j, o~ y j, 
and 



are the the standard deviations of the x, y, 



and y j3 vectors over the full range in i. In this case the 
standard deviation does not represent the uncertainty on 
the measured positions but instead reflects the relative 
range in position spanned by the points in our vector. By 
using this parameter to determine our smoothing width, 
we effectively implement an adaptive smoothing scheme 
with a smaller spatial resolution in regions with dense 
sampling and a larger resolution in regions where the 
pixel is more sparsely sampled. The i index sums over 
the nearest 50 neighbors with distance dij defined as 
follows: 



dij = (a 



Vi ~ Vj 



+ V A - ( 3 ) 



where b is an empirically determined scaling constant 
that allows more weight to be placed on either the x or 
y separation. Both b and the total number of points n 
used in Eq. [3] are fixed at the values that produce the 
lowest standard deviation in the residuals from the best- 
fit solution. In this case we obtain the best results when 
b = 0.8 in the 3.6 /jm band, while we prefer b = 1.0 in 
the 4.5 /im band. We find that increasing the number of 
nearest neighbors used in the weighted average from 50 
to 400 reduces the relative scatter in the final unbinned 
residuals by less than 1% and increases the computation 
time by a factor of six, making our Markov chain fits 
computationally intractable. 

We calculate the intrapixel sensitivity correction using 
the residuals after the best fit functions for the phase 
curve (£ )2.4I) . transit and secondary eclipses f ^2.3p . and 
stellar activity f £12.51) have been removed in order to en- 
sure that we do not inadvertently remove part of these 
astrophysical signals. We re-calculate this correction at 
each step in our fits, allowing the pixel map to trade off 
against changes in our nominal astrophysical functions. 
We find that this iterative approach results in a unique 
best-fit solution despite the large amplitudes of the in- 
trapixel sensitivity variations, as improved estimates for 
the astrophysical parameters result in a more accurate 
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pixel map, which in turn results in a lower scatter in the 
residuals after both the astrophysical signals and the in- 
trapixel sensitivity variations have been removed. It is 
worth noting, however, that the success of this approach 
depends on our ability to generate accurate models for 
the astrophysical signals in our data. 

2.3. Transit and Eclipse Fits 

We calculate our eclipse curve using the equations 
from IMandel fc Agoll (|2002ft assuming a circular orbit 
for the planet, as ind i cated by both radial velocity data 
(IBouchv et all 120051: iTriaud et all 120091: IBoisse et af 
2009) and secondary ec lipse photometry (jAgol et al 
2010 : Ide Wit et all 12012ft. We f ix the planet's period 



to the value from iKnutson et al.l (|200 9a). and tak e our 
limb-d arkening coefficients in each bandpass from ISinri 
(20100 assuming an ATLAS stellar atmosphere model 
with T eff = 5100 K, log(p) = 5.0, and [Fe/H]=0 
(jBouchv et al J 12005ft . Our transit light curve has four 
free parameters, including the planet-star radius ratio 
R P I R+, the orbital inclination i, the ratio of the orbital 
semi-major axis to the stellar radius a/i?*, and the tran- 
sit time. For the secondary eclipse we appro ximate the 
plane t as a uniform disk (for more on this see lAgol et al.1 
12010ft where the depth is measured relative to the av- 
erage of the phase curve function over the duration of 
the eclipse. We account for the change in the planet's 
brightness during the secondary eclipse by rescaling the 
amplitudes of ingress and egress to match the values of 
the phase curve at the start and end of the eclipse, re- 
spectively. Our secondary eclipse light curves use the 
same values of R p /R+, i, and a/R* as the transit, and 
we allow the depths and times of each secondary eclipse 
to vary individually. This gives us a total of eight free 
parameters in our fits to describe the transit and two 
secondary eclipses. 

2.4. Phase Curve Functions 

We calculate our phase curve functions as the sum of 
a series of sine and cosine terms with periods equal to 
the planet's orbital period and i ts even harmonics as 
described in iCowan fc Agoll (|2008ft . This approach pro- 
duces phase curves with shapes that are indistinguishable 
from those of the "oran ge-slic e" models we used previ- 
ously in IKnutson et al.l (|2007l l2009aft . but with the ad- 
vantage that they result in more physically motivated 
smoothly varying brightness distributions instead of dis- 
crete segments of constant brightness on the surface of 
the planet. We tried fits including sine and cosine terms 
with periods equal to the planet orbit P, as well as 
higher-order harmonics of P/2, P/4, and P/8. As a 
check we also tried fitting the o dd harmonics P/3 an d 
P/6 to the data; as discussed in ICowan fc Ag ol (2008), 
these harmonics cannot be produced by a static plane- 
tary brightness map, and their amplitudes therefore pro- 
vide a useful test for whether or not we are over-fitting 
the data. We find that the odd harmonic terms have fit- 
ted amplitudes comparable to those of the P/4 and P/8 
terms, and therefore elect to limit our fits to the P and 
P/2 harmonics: 

13 Available at www.astro.ex.ac.uk/people/sing 
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Figure 3. Visible-light flux variations measured for HD 189733 in 
the averaged Stromgren 6 and y filters (black circles) with the APT 
spanning the epoch of our 3.6 (im phase curve observation (vertical 
blue lines). The solid red line indicates the best-fit spot solution 
with lc uncertainties delineated by the shaded pink region, which 
we use to extrapolate the likely stellar flux variations during our 
phase curve observation. 
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Figure 4. Visible-light flux variations measured for HD 189733 in 
the averaged Stromgren b and y filters (black circles) with the APT 
spanning the epoch of our 4.5 fxra phase curve observation (vertical 
blue lines). The solid red line indicates the best-fit spot solution, 
with lc uncertainties delineated by the shaded pink region, which 
we use to extrapolate the likely stellar flux variations during our 
phase curve observation. 



F{\) = l + CiCOs(27rf/P) 
+C3 cos(47rf/P) 



c 2 sin(27rf/P) 
c 4 sin(47rf/P) (4) 



where c\ — C4 are free parameters in the fits and t is 
the time in days from the predicted center of transit. 

2.5. Stellar Variability 

HD 189733b is an active K0 star with a Ca II H fc K 
emission line strength of lo g(P^ T< -)= —4.5 (|Bouchv et all 
120051 : IKnutson et al.ll2010ft . which has been observed to 
vary by as much as ±1.5% at visible wavelengths with 
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Bin Size 

Figure 5. Standard deviation of residuals vs. bin size for the 
3.6 fim (top panel) and 4.5 fira (bottom panel) data after correct- 
ing for intrapixel sensitivity variations and dividing out the best-fit 
solutions for the transit, secondary eclipses, phase variations, and 
stellar activity. Bin sizes of 10 3 and 10 correspond to time inter- 
vals of approximately 2 minutes and 20 minutes, respectively. 

a rotation per iod of 11.95 days (iHenry fe Winnl 120081 : 
iKnutson et aLll2009al lLanza et al.l 120111 ). Although the 
amplitudes of these variations are reduced in the in- 
frared, they are still comparable to the planet's phase 
curve and must be account ed for in our estimat es of 
the phase curve amplitude (jKnutson et al.l l2009af) and 
wavelength-dependent transit depth (iPont et al.1 120081 : 
IDesert et al.l [20091 [20lli ISing et al.l 1201 ID . We char- 
acterize the nature of these trends using ground-based 
monitoring data for HD 189733b obtained with the Ten- 
nessee State University T10 0.8 m automated photomet- 
ric telescope (APT) at Fairborn Observatory as part of 
a long-term monit oring program extend over a period of 
six years m total (|Henrv fc Winnl 120081) . The data were 
obtained in the Stromgren b and y filters, nodding be- 
tween HD 189733 and three comparison stars of compa- 
rable o r greater brightness as described in lHenrv fc Winnl 
(2008). We present the averaged photometry spanning 
the epochs of the 3.6 and 4.5 /jm Spitzer data in Fig. |3] 
andgl 

Because HD 189733 was not observable from the 
ground during the dates of the Spitzer observations, we 
must extrapolate from ground-based measurements ob- 
tained several weeks prior and following our observa- 
tions in order to infer the behavior of the star during 
this period. We fit the data shown in Fig. [3] and 0] 
with a spo t model for th e star following the method de- 
scribed in lAigrain et al.l ([201 ll ). with uncertainties cal- 
culated using a Gaussian processes method described 



in iGibson et al.l (|2011bl ). Based on these models, we 
find that the star varied in brightness by approximately 
0.4%- 1.8% (3.6 /im) and 0.3%- 1.4% (4.5 /on) in the av- 
eraged Stromgren b and y filters during our phase curve 
observations, with a local minimum in the 3.6 /im data 
and a linearly decreasing flux in the 4.5 /im data. If we 
rescale these variations to the wavelengths of our Spitzer 
obse rvations assumin g a spot temperature of 4250 K 
(e.g. jPont et~aL|[2008l) . we find corresponding amplitudes 
of 0.1% - 0.4% at 3.6 /mi and 0.1% - 0.5% at 4.5 /on. 
We account for the effects of the spots in our data by 
including a quadratic function of time in our 3.6 /im fit 
and a linear function of time in our 4.5 /im fit: 

F(t) = dit + d 2 t 2 (5) 

where t is defined as the time from the predicted center 
of transit, d\ and d 2 are free parameters in the 3.6 /jm fit, 
and we set d 2 equal to zero in the 4.5 /im fit. As a test 
of our ability to distinguish between spot variability and 
the planet's phase curve in these bands, we inserted lin- 
ear functions of time with slopes ranging between ±0.5% 
in both data sets. We find that we are able to retrieve the 
correct slope for our inserted trend in all cases, although 
our nominal phase curve solutions displayed some degen- 
eracies with these trends for amplitudes larger than 0.3%. 
Including a quadratic function of time in our fits to the 
the 4.5 /jm data produces a negligible improvement in 
the Bayesian Information Criterion (BIC; 1,651,608 for 
the linear fit vs. 1.651.584 for the quadratic fitful. As an 
additional test, we generated 100 random samples from 
the Gaussian process model trained on the visible-light 
photometry. We then fit the predicted stellar trend from 
each sample during the epoch of our 4.5 /im phase curve 
observation with both a linear and a quadratic function 
of time. Wc find that the addition of a quadratic term 
produces a significant (ABIC > 5) improvement less 
than 30% of the time, in good agreement with our pre- 
vious test. We therefore limit our 4.5 /tm fits to a linear 
function of time in the analysis described below and find 
that the star decreased in brightness by 0.180%±0.022% 
over the course of our observations (see Table[T]), in good 
agreement with our prediction of 0.1% — 0.5% from the 
visible-light photometry. 

In contrast to this result, we find that at 3.6 /tm the 
addition of a quadratic term to our stellar spot func- 
tion results in a significantly improved BIC (1,604,852 
for the linear fit vs. 1,604,538 for the quadratic fit). 
This is consistent with the results from our sampling 
of the Gaussian process model, which indicate that the 
quadratic term provides a better description for the pre- 
dicted trend during our 3.6 /im phase curve observations 
in more than 70% of cases. However, we find that the 
best-fit quadratic solutions from our fits to the visible- 
light and 3.6 /im photometry are mutually inconsistent; 
we strongly prefer a negative quadratic term and a local 
flux maximum in our 3.6 /im fits, whereas the visible-light 
photometry indicates a positive quadratic term and a lo- 
cal flux minimum in 83% of our solutions. Upon closer 
examination we noticed that our best-fit quadratic func- 
tion of time from the 3.6 /tm photometry produced a 

14 Defined as BIC= x 2 + feln(n), where k is the degrees of 
freedom and n is th e total number of points in the fit (e.g. 
IStevenson et al.ll20Tch . 
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Figure 6. Final 3.6 fim photometry after correcting for in- 
trapixel sensitivity variations and stellar activity (black filled cir- 
cles), binned into four minute intervals. The best-fit phase, transit, 
and eclipse curves are overplotted as a red line. The lower panel 
shows the same data as the upper panel, but with an expanded y 
axis for a better view of the phase curve. 

local maximum just before the secondary eclipse, around 
the same time where we would expect to see a maxi- 
mum in our best-fit phase curve. We therefore conclude 
that the quadratic fit is degenerate with our phase curve 
function in this bandpass and opt to limit the quadratic 
coefficient in our subsequent fits to positive values (i.e., 
a linear function or a local flux minimum) in order to 
avoid this degeneracy and ensure good agreement with 
our visible-light spot models. Our best-fit solution in this 
band is very close to linear and indicates that the star 
increased in brightness by 0.262% ± 0.016% during our 
observations (see Table [TJ in good agreement with our 
prediction of 0.1% — 0.4% from the visible-light photom- 
etry. 

3. RESULTS 

We fit the trimmed data and calculate uncertainties 
simultaneously for the transit, secondary eclipses, phase 
curve, stellar activity, and intrapixel sensitivity cor- 
rections using a Markov Chain Monte Carlo (MCMC) 
method (see, for exa mple Ford 2005t iWinn et al.l l2007b) 
with a total of 10 5 steps and either 16 (3.6 fim) or 
15 (4.5 /im) free parameters. These parameters in- 
clude: the four phase function coefficients c\ — C4, a/R±, 
i, Rp/R+, transit time, two secondary eclipse depths, 
two secondary eclipse times, either a linear (4.5 /jm) or 
quadratic (3.6 fim) function of time to account for stel- 
lar variability, and two noise parameters discussed be- 
low. We plot the normalized time series after the best-fit 
intrapixel sensitivity variations and stellar trends have 
been removed in Fig. |6]and[7] 

We find that there is still some remaining time- 
correlated noise in the residuals from our best-fit solu- 
tion (Fig. [5]); we account for this time correlation by 
im plementing a wavelet-b ased MCMC fit as described 
in iCarter &: Winnl (|2009f ) and compare the results of 
this fit to t he st andard y 2 -based m ethods described in 
iFordl (|2005D and IWinn etaTI (|2007bl) . For the standard 
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Figure 7. Final 4.5 fim photometry after correcting for intrapixel 
sensitivity variations and stellar activity (black filled circles); see 
Fig. [6] for a complete description. 

MCMC fit, which assumes that the uncertainties on in- 
dividual points are Gaussian and time-independent, we 
set the per-point uncertainty equal to the value needed to 
produce a reduced x 2 equal to one for the best-fit solution 
(0.333% at 3.6 fim and 0.466% at 4.5 fim). In our wavelet 
MCMC we maximize the likelihood function instead of 
minimizing % 2 , which allows us to fit for the white (a w ) 
and red (ay) noise contributions to the per-point uncer- 
tainties. The wavelet transform routine requires our data 
to have a length equal to a power of two, which we achieve 
by subdividing each phase curve into thirteen segments 
of equal length and zero-padding by either 6% (3.6 /im) 
or 3% (4.5 /im). We prefer zero-padding over trimming 
as it allows us to include all available data in our fits, 
and we find that this approach has a negligible effect on 
our best-fit noise parameters. 

We initialize each chain at a position determined by 
randomly perturbing the best-fit parameters from a 
Levenberg-Marquardt minimization. After running the 
chain, we search for the point where the \ 2 value first 
falls below the median of all the % 2 values in the chain 
(i.e. where the code had first found the optimal fit), and 
discard all steps up to that point. For the wavelet MCMC 
we perform a similar trim where the likelihood first rises 
above the median value. We set our best-fit parameters 
equal to the median of each distribution and calculate 
the corresponding uncertainties as the symmetric range 
about the median containing 68% of the points in the 
distribution. The distribution of values was very close 
to symmetric for all parameters, and there did not ap- 
pear to be any detectible correlations between variables 
aside from i and a/i?*. Our wavelet analysis produces 
modestly higher uncertainties for most parameters, with 
the greatest increases in the planet-star radius ratio and 
secondary eclipse depth errors, which were a factor of 
two larger than in the standard Gaussian \ 2 fits. We list 
the best-fit parameters and corresponding wavelet-based 
uncertainties in Table [T] 

For both fits the standard deviation of our best-fit 
residuals is a factor of 1.12 higher than the predicted 
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Table 1 

Global Fit Parameters 



Parameter 



3.6 /im 



4.5 /im 



Transit Parameters 
Rpl -R* a 

T c (BJD) b 

Secondary Eclipse Parameters 

1 st Eclipse Depth 

T c (BJD) b 

2 nd Eclipse Depth 

T c (BJD) b 

Phase Curve Parameters (Eq. 

ci (cos(27rt/P) term) 

C2 (sin(27rt/P) term) 

C3 (cos(47rt/P) term) 

C4 (sin(4-7rt/P) term) 

Amplitude 

Minimum Flux Offset c [h] 
Maximum Flux Offset [h] 

Stellar Flux Variation (Eq. 
Linear term di 
Quadratic term d^ 

Noise Parameters 

o~ w 
^2e 



0.15511 ±0.00020 
85.671 ±0.034 
8.858 ±0.031 
2455559.554550 ± 0.000035 



0.1440% ±0.0053% 
2455558.44753 ± 0.00081 

0.1500% ±0.0061% 
2455560.66515 ± 0.00046 



-0.0479% ± 0.0032% 
-0.0389% ± 0.0028% 
0.0025% ± 0.0021% 
-0.0020% ± 0.0014% 
0.1242% ±0.0061% 
-6.43 ± 0.82 
-5.29 ±0.59 



0.000932 ± 0.000057 
< 0.0000041 



0.0032608 ± 0.0000027 
0.000581 ± 0.000043 



0.15580 ±0.00019 
85.735 ±0.036 
8.902 ±0.032 
2455189.052491 ± 0.000032 



0.1889% ±0.0056% 
2455187.94447 ± 0.00040 

0.1696% ±0.0053% 
2455190.16158 ±0.00048 



-0.0468% ± 0.0044% 
0.0122% ± 0.0027% 
-0.0011% ±0.0019% 
-0.0020% ± 0.0022% 
0.0982 ± 0.0089% 
-1.37 ± 1.00 
-2.98 ±0.82 



-0.000652 ± 0.000079 
(fixed) 

0.0046009 ± 0.0000037 
0.000283 ± 0.000037 



a These values do not include a correction for stellar variability; sec ^4.21 for the corrected values. 
b We list all times in BJD_UTC for consistency with other studies; to convert to BJD_TT, simply add 
66.184 s. 

c Offsets are measured relative to the center of transit time for the minimum flux and the center 
of secondary eclipse time for the maximum flux. A negative offset indicates that the maximum or 
minimum occurs earlier than predicted, corresponding to an eastward offset in the location of the hot 
or cold region in the planet's atmosphere. 

d This parameter is limited to values greater than or equal to zero in our fits; sec SI2.5I for more 
information. We find that our best-fit distribution peaks at zero and we therefore report the la upper 
limit here. 

e Unlike a w , the red noise parameter is not equal to the variance of the red noise component of the 
time scries; sec ICarter fc Winnl 120091) for a more detailed explanation. 



photon noise limit at 3.6 /im, and 1.13 times higher than 
this limit at 4.5 /im. If we calculate the standard de- 
viation of the red noise vector in our best-fit wavelet 
solution, we find that the red noise contributes 0.0162% 
(3.6 /im) and 0.0017% (4.5 /im) of the total scatter in the 
unbinned residuals. We also estimate the red noise on 
time scales relevant for the eclipses (i.e., approximately 
one hour) by comparing the standard deviation of our 
best-fit residuals to the theoretical Poisson limit for a 
bin size of 27,000 points in Fig. [5j We find that red 
noise contributes a scatter of approximately 0.003% in 
relative flux in both bandpasses, corresponding to a fac- 
tor of 2.5 (3.6 /im) or 2.1 (4.5 /im) increase in the total 
noise for our one hour bins. This is consistent with the 
results of our wavelet analysis, which finds that the un- 
certainties in the measured secondary eclipse depths in 
both channels are approximately double those predicted 
assuming Gaussian noise in the standard MCMC fits. 

4. DISCUSSION 

We present the results of our transit and secondary 
eclipse fits in the sections below, focusing on the mea- 
sured depths and times of eclipse. We find an average 
orbital inclination of 85.701° ± 0.025°, a/R„ = 8.879 ± 
0.022, and an impact parameter b = 0.6656±0.0042, con- 



sistent wit h prev io us estimate s from iPont et all (12 0071 , 
iSing et aU (120091) . lAgol et all (|2010l ). and IDesert et all 
(2009, 2011a). We note that our errors are smaller 
than those reported by IDesert et al.l (|2011a| ) for their 
3.6 /im transit observation; this is likely due to the avail- 
ability of a longer baseline and correspondingly tighter 
constraints on the correction for the intrapixel sensitivity 
variations in our observations. 

4.1. Orbital Ephemeris 

We show the best-fit transit and secondary eclipse so- 
lutions in Fig. [8j we see no evidence for excess scat- 
ter in the transit residuals, suggesting that our assump- 
tions for the shape of the planet (i.e., oblate vs. circular) 
and limb-darkening coefficients are correct to within our 
uncertainties. We calculate an updated ephemeris for 
the system using the transit times given in Table [J and 
previously published values from Spitzer dKnutson et al.l 
l2009at IDesert et all 120091 12011a!: lAgol et al.l 120101) and 



Hubble (jPont et al.ll2007l ). We define this ephemeris as: 



T c (n) = T c (0)+n x P 



(6) 



where T c is the predicted transit time as a function of 
the number of transits elapsed since T c (0) and P is the 
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Figure 8. Best-fit transit and secondary eclipse light curves for the 3.6 (im (a— f) and 4.5 /an (g— 1) bands (black filled circles). The data 
have been binned in four minute intervals. The first secondary eclipses are shown on the left (a and g), the second secondary eclipses are 
shown in the middle (b and h), and the two transits are shown on the right (c and i). Residuals from the best-fit solution (black line, a— c 
and g— i) are shown immediately beneath each event in panels d— f and j— 1. 



orbital period. Our new observations extend the pre- 
vious baseline by more than two years, and as a result 
we derive a more precise estimate of the planet's orbital 
ephemeris: T c (0) = 2455559.554587 ± 0.000025 BJD and 
P = 2.218575143 ± 0.000000046 days (±4 ms), where 
T c (0) and P are determined from a linear fit to our mea- 
sured transit times and the previously published values 
listed above. As shown in Fig. [9j our new transit times 
are consistent with a constant ephemeris in all cases. 

We fit a second, indepen dent ephemeris to t he sec - 
ondary eclipse times from iCharbonneau et "all ([2008) , 
iKnutson et ail (|2009aft . and lAgol et all (|2010D . as well 
as this paper, and find a best-fit value of: To = 
2455560.66451 ± 0.00025 BJD and P = 2.218574560 ± 
0.000000454 days. This period differs from the best-fit 
transit period by 1.3er. 

We find that the secondary eclipse times occur sys- 
tematically later than predicted, with an average de- 



lay of 49 ± 9 s after accounting for the 31 s. light 
travel time delay (jAgol et al.ll2010f ). The additional de- 
lay is consistent with the time offset predicted for a 
planet with a non-uniform surface brightness, which al- 
ters the shape of ingress and eg ress and the correspo nding 
best-fit ecl i pse tim e estimate (| Williams et al.ll2006[ ). In 
lAgol et all j|010) we found that the eastward offset of 
the hot spot on the planet's day side causes the best-fit 
8 /zm eclipse time to shift approximately 30 s later than 
the true valu e . Th is res ult was recently co nfirmed by 
iMaieau et all (|2012fl and (|de Wit et al.ll2012fl . who used 
a variety of methods to create 2D maps of HD 189733b's 
8 /xm dayside brightness distribution based on the shape 
of its 8 /xm phase curve and secondary eclipse ingress 
and egress. If we examine our new secondary eclipse 
times separately we find average offsets of 99 ± 35 s at 
3.6 fim and 22 ± 27 s at 4.5 /mi after accounting for the 
light travel time. Although these estimates are consis- 
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Figure 9. Observed minus calculated transit (top panel) and sec- 
ondary eclipse (bottom panel) times from the published literature. 
The two most recent transit measurements and four most recent 
secondary eclipse observations are derived from the observations 
presented in this paper; the timing offsets in earlier transits are 
likely the result of occulted star spots. We calculate predicted sec- 
ondary eclipse times assuming an offset of half an o rbit from the 
transit center and use the ephemeris presented in EI4.1I Dashed 
lines indicate the ltr uncertainties in the predicted transit and sec- 
ondary eclipse times. 

tent with the 8 /jm offsets, they are not precise enough 
to provide independent constraints on the location of the 
dayside hot spot in the 3.6 and 4.5 /an bands. 

4.2. Transit and Secondary Eclipse Depths 

We obtain planet-star radius ratios of 0.15511±0. 00020 
at 3.6 jam and 0.15580±0. 00019 at 4.5 /am, which we com- 
pare to previously published values of 0.154 5 ± 0.0003 at 
3.6 /zm and 0.1557 ± 0.0003 at 4.5 /zm froml Desert et all 
( 2"509l) and O.^ee ^^ at 3.6 /zm from Desert eTaTI 
( 2011af) . Although our results are consistent with these 
previous measurements at the 1.7a (3.6 /zm) and 0.3cr 
(4.5 /im) level, we note that variations in the stellar flux 
between epochs can affect the measured transit depths, 
and must be accounted for if we wish to compare these 
two sets of measurements. In this case we are most in- 
terested in the relative difference between our measured 
3.6 and 4.5 /zm radius ratios, which can be used to search 
for the presence of absor ption feature s in HP 189733b's 
transmission spectrum. IDesert et al.l ()2011al ) find that 
their 4.5 /im radius ratio is 0.0012 ± 0.0004 smaller than 
at 3.6 /zm and suggest that this may be consistent with 
absorption from carbon monoxide, although they do not 
correct explicitly for changes in the stellar flux over the 
approximately one year interval between their two mea- 
surements. 



We account for the effect of varying stellar flux in our 
data using the ground-based monitoring data described 
in g2T5j We find that the star is 0.9%±0.6% fainter in the 
Stromgren (b + y)/2 band during our 3.6 /zm transit as 
compared to our 4.5 /zm transit. If we assume a spot tem- 
perature of 4250 K as estimated bv lPont e t al. (2008|), we 
find that the corresponding flux change in the 3.6 and 
4.5 /tm bands should be approximately 28% of the vari- 
ation in the Stromgren (b + y)/2 band, or 0.3%±0.2%. 
Correcting the 3.6 /zm transit depth to the equivalent 
stellar flux at the time of the 4.5 /zm transit observa- 
tion yields a planet-star radius ratio of 0.15465±0. 00037, 
0.001 15±0.00042 (2.7a) smaller than in the 4.5 /tm band. 
This result is consistent with the planet-st ar radius ra- 
tio dif ference of 0.0012 ± 0.0004 reported by IDesert et~aTl 
(2009), who also find a deeper transit at 4.5 /zm as com- 
pared to 3.6 /zm. We note that although our measure- 
ments of the planet-star radius ratios are more precise 
than those reported in this paper, we are ultimately lim- 
ited by the uncertainties in our estimates of the stellar 
flux and its corresponding effect on the measured transit 
depths at the epochs of our two measurements. 

We also present new estimates of the secondary eclipse 
depth in the 3.6 and 4.5 /zm bands. We see no con- 
vincing evidence for variability in the planet's flux over 
the planet's orbit; our two 3.6 /zm depths differ by 0.8cx, 
and our 4.5 /zm depths differ by 2.5rj. We estimate an 
average secondary eclipse depth of 0.1466% ± 0.0040% 
at 3.6 /tm and 0.1787% ± 0.0038% at 4.5 /tm. We uti- 
lize a PHOENIX stellar atmosphere model for the star 
and calculate corresponding brightness temperatures of 
1328 ± 11 K and 1192 ± 9 K, respectively. These val- 
ues are shallower than our previous estimates of 0.256 ± 
0.014% at 3.6 zzm and 0. 214 ± 0.020% at 4.5 /zm from 
iCharbonneau et all (|2008l ). although the 4.5 /zm depths 
are consistent at the 1.7tr level. We note that the er- 
ror estimates for our new secondary eclipse depths are a 
factor of 4 — 5 smaller than these previous values; this 
is mainly due to the fact that the older observations cy- 
cled between all four Spitzer IRAC bands in order to ob- 
tain multi-wavelength coverage during a single secondary 
eclipse event, leading to an effective cadence one- fifth 
that of our new staring-mode observations. Although 
our new 3.6 /zm secondary eclipse depth differs from our 
older estimate by 7.5ct, the older observation had much 
larger intrapixel sensitivity variations that were partially 
degenerate with the secondary eclipse depth in this band. 
Because our error estimates in the 2008 paper were based 
on conventional MCMC methods, which assume that 
the noise is Gaussian, it is likely that the uncertainties 
on the older eclipse depth were significantly underesti- 
mated. Although we considered re-analyzing these older 
data with our new methods, we concluded that the lower 
cadence and larger pointing jitter in these observations 
meant that we were unlikely to get useful constraints on 
the secondary eclipse depth as compared to our new ob- 
servations. 

4.2.1. Comparison to ID Atmosphere Models 

We compare our new secondary eclipse depths to the 
predictions of one-dimensional radiative transfer mod- 
els for HP 18 9 733b's dayside emission as described in 
iBurrows et al.l (|2008l ). These models assume a solar 
metallicity atmosphere with the chemistry in local ther- 
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mal equilibrium (LTE), an interpolated ATLAS model 
for the stellar spectrum, and allow for the parameter- 
ized redistribution of heat from the dayside to the night 
side using a heat sink between 0.01 — 0.1 bars. They 
also include an additional grey absorber at low pres- 
sures, which is required to reproduce the high-altitude 
temperature inversions observed in the atmospheres of 
a subset of hot Ju p iters such as HP 2 9458b (e.g.. 
Burrows et all 120071: IKnutson et al.l 120081: iSwain et alj 
20091 : IMadhusudhan fc Seagerll2010D . Both the fractional 
amount of the incident flux P n redistributed to the night 
side and the opacity of the high-altitude absorber n e 
in cm -1 g are left as free parameters and optimized 
to provide the best match to the updated secondary 
eclipse depths at 3.6 and 4.5 fim, together with previ- 
ously pubjished_5 ; 8 J _8 ; 0, 16, an d 24 urn eclipse depths 



from iCharbonneau et"afl (j2008j ) and 
and th e 5—12 /im IRS spectrum from iGrillmair et alJ 
( 2008). We obtain a reasonably good match to all mea- 
surements except the 3.6 /j,m dayside using a model with 
P n = 0.15 (i.e., 15% redistribution, where 50% redis- 
tribution corresponds to a perfectly well-mixed atmo- 
sphere) and K e — 0.035 cm -1 g, as shown in Fig. [TUJ 
These values are identical to those of the nominal model 
presented in IGrillmair et all (|2008l ). For this model, our 
measured dayside fluxes differ by [— 12.1<7, — 0.9<r, +1.6cr, 
—O.lcr, +3.6(7, +2.9er] from their predicted values in the 
[3.6, 4.5, 5.8, 8.0, 16, 24] ^m bands, respectively. 

Given that our models include an estimate of the rel- 
ative day-night recirculation efficiency, it is worth con- 
sidering whether or not a nightside model with the same 
assumptions as our day side model can accurately re- 
produce the measured night side fluxes at 3.6, 4.5, 8.0, 
and 24 ^m. We take the l onger- w avelength flux es- 
timates from IKnutson et alJ ([2009a) and estimate the 
3.6 and 4.5 /Ltm night side fluxes as the measured sec- 
ondary eclipse depths minus the phase curve ampli- 
tude in each band; this gives planet-star flux ratios of 
0.0224% ± 0.0073% and 0.0805% ± 0.0097% at 3.6 and 
4.5 /im, respectively. If we use a PHOENIX stellar at- 
mosphere model for the star, these fluxes correspond to 
brightness temperatures of 825 ± 54 K and 928 ± 33 K, 
respectively. We compare these fluxes to the predictions 
of a solar metallicity LTE model with P n — 0.15; the 
model does not include any additional high-altitude ab- 
sorber, but this is irrelevant for the night side as there is 
no incident flux at the top of the atmosphere. We find 
that this default night side model provides a good match 
for the measured 3.6 and 4.5 /jm nightside fluxes (—1.9(7 
and +2.3(7 difference, respectively), although it underes- 
timates the 24 fim flux by 7.1a. We conclude that statis- 
tical surveys seeking to estimate the relative amount of 
day-night circulation based on fits to th e dayside emis- 
sion a lone, such as those presented in ICowan fc Agoll 
(|2011l ). should obtain reasonable estimates for this quan- 
tity despite the limited data available and the relative 
simplicity of the models used. 



4.2.2. A Note on Stellar Atmosphere Models 

There are several potential sources of stellar atmo- 
sphere models for HD 189733 in the literature, includ- 
ing both interpolated and star-specific ATLAS models 




Wavelength (micron) 



Figure 10. Upper panel: measured broadband secondary eclipse 
depths (filled da rk blue circles) at 3.6 and 4.5 fim (this paper), 5.6, 
16, a nd 24 (im {Charbonncau et al. 2008), and 8 /mi lAgoI et alj 
2010), and an IRS spectrum (open light blue circles, Grillmair et al. 
2008) compared to a ID atmosphere model fo r HD 189733b's day 
side (black line) following Burrows et al. (2008) with parameterized 
recirculation P n = 0.15 at depth and a high-altitude absorber with 
an opacity equal to K e = 0.035 cm -1 . Black filled squares show the 
predicted planet-star flux ratios in the same bands as the Spitzer 
observations. The lower panel shows the measured minimum fluxes 
for HD 189733 (filled dark blue circles) compared to the same ID 
atmosphere model for the planet's night side; the temperature in 
this region is set by the amount of energy transported from the 
day side. 



(!Kuruczlll979l. 11994 Effl OSfl and interpolated PHOENIX 
models (jHauschildt et al.1119991) . We compare the stellar 
spectra obtained from these three sources, and find that 
the star-specific ATLAS model over-predicts the strength 
of the CO absorption feature between 4 — 5 fim as com- 
pared to our interpolated ATLAS and PHOENIX spectra. 
If HD 189733's spectrum is best-described by this star- 
specific model it would shift the predicted planet-star 
flux ratio in the 4.5 /xm Spitzer band to larger values, 
at a level comparable to or larger than our uncertain- 
ties for this measurement. For now we proceed using the 
interpolated grid models, which we find give consistent 
predictions regardless of which grid they are taken from. 

4.3. Phase Curve Fits 

There are three basic quantities which we obtain from 
our phase curve observations: the amplitude of the ob- 
served variation, the location of the flux maximum, and 
the location of the flux minimum. Taken together, these 
quantities can be used to estimate the day-night bright- 
ness contrast as well as the positions of hot and cold 
regions in the atmosphere. Because the opacities (and 
hence the pressure from which most of the detected flux is 
emitted) may vary from one band to the next, we do not 
necessarily expect to see the same features in all band- 
passes. This picture is further complicated by the fact 

15 Available a.t \http : / /kurucz. cfq.ha rvard.edu/ 
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that the wavelength-dependent opacity of the planet's 
atmosphere may vary as a function of longitude. Our 
first task, therefore, is to consider whether or not our 
ensemble of phase curves all display the same qualitative 
features. 

Focusing on the most basic feature first, we compare 
the amplitude of the phase variation in our four bands. 
We find that the planet's minimum flux is 15%±5% of the 
maximum flux at 3.6 /im, and 45% ± 5% at 4.5 /tm. We 
can then convert this flux contrast to a brightness tem- 
perature contrast of 503 ±21 K at 3.6 /mi and 264 ±24 K 
at 4.5 /im. The contrast in the 3.6 /tm band is 6 — 7a 
larger than in the 4.5 /im ba nd or the correspondin g value 
of 188±48 K at 24 /im from lKnutson et al.l p009al) . If we 
assume that the composition of the planet's atmosphere 
is uniform (e.g., the opacity is constant), this would sug- 
gest that at longer wavelengths we are seeing deeper 
into the planet's atmosphere, where t emperatures are 
expected to be more homo genous (e.g... iShowman et al.l 
[20091: iBurrows et all [2010h . This differs from the pre- 
dictions of ID solar metallicity at mosphere models as- 
suming equilibrium c hemistry (e.g.. iFbrtnev et al.l 120081 : 
IBurrows et"all 120081 : Knutson et al l l2009af l . which find 
that HD 189733b's 3.6 /im band should have the low- 
est average opacity and, therefore, probe deepest in the 
planet's atmosphere. 

We next compare the locations of the flux maxima in 
each band. If the planet's atmosphere is in strict radia- 
tive equilibrium (i.e., no recirculation) we would expect 
the phase curve to reach its maximum at the center of 
the secondary eclipse, when we are viewing the substcl- 
lar point face-on. If, on the other hand, the atmospheric 
circulation is dominated by a strong west-to-east (i.e., 
superrotating) equatorial band of wind, we would expect 
to see the phase curve reach a maximum prior to the 
secondary eclipse. We find that the maximum flux oc- 
curs 5.29 ± 0.59 hours before the center of the secondary 
eclipse at 3.6 /tm and 2.98 ±0.82 hours before the eclipse 
center at 4.5 /tm. These two offsets differ from zero by 
9.0 and 3.6<7, respectively, and from each other by 2.3cr. 

We next compare the locations of the minima in both 
light curves, which should be offset in the same direction 
as the maxima in the limit of a single equatorial band of 
wind. We find that the minimum flux occurs 6.43 ± 0.82 
hours before the transit center at 3.6 /im and 1.37 ± 1.00 
hours before the transit center at 4.5 /tm. These two 
measurements differ from zero by 7.8 and lAa respec- 
tively and are consistent with the presence of an east- 
ward jet on the night side. If this same eastward offset 
was present in our 8.0 and 24 /tm phase curves it would 
have gone undetected, as our observations began shortly 
before the transit and spanned just half an orbit in both 
cases. General circulation models ( e.g.. IShowman et al.1 



120091: IBurrows et al.ll2010HHeng et al.ll2011h predict that 
in cases where a super-rotating eastward jet dominates 
heat transport the cold region on the planet's night side 
should display a larger shift than the hot region on the 
dayside; this results in maximum and minimum fluxes in 
the integrated light curves that are separated by slightly 
less than 180°. We find that our data are consistent with 
a 180° separation in both bands, but this is not surpris- 
ing as our constraints are weak compared to the size of 
the predicted deviations. 



4.3.1. The Location of the 8 fim Phase Curve Minimum 

Spitzer data taken in the 8 /tm band display a ramp- 
like behavior where higher-illumination pixels converge 
to a constant value within the first hour of observa- 
tions while lower-illumination pixels show a continually 
increasing linear trend on the time scales of interest here. 
This effect is generally attributed to cha rge-trapping in 
the array an d is discussed in deta il in Knutson et al.l 
([20071 l2009ch and lAgol et all (pOlOl ). among others. We 
must correct for t his effect in order to re trieve the phase 
curve presented in lKnutson et al.l (|2007| ). which displays 
a flux minimum several hours after the transit (i.e., 
shifted in the opposite direction of the minima in our 
new 3.6 and 4.5 /im phase curves). Although we initially 
considered this minimum to be reliably detected despite 
the uncertainties in o ur ramp cor r ection , an independent 
analysis presented in lAgol et all (|2010| ) used a different 
method to correct the ramp and found much larger un- 
certainties at early times. 

In this section we investigate the nature of the detector 
ramp in more detail, focusing on the question of whether 
or not the flux minimum in the 8 /im data might be ei- 
ther a direct result of the ramp or an artifact created by 
our ramp correction. We utilize 8 /tm preflash data ob- 
tained on U T 2009 Apr 23 as pa rt of a program to observe 
CoRoT-7b (|Fressin et al.ll2012fl . which targeted a region 
in the Orion Nebular Cloud with bright, diffuse emission 
located at J2000 05:35:16 -05:23:24. The purpose of these 
observations was to mitigate the ramp effect by filling 
the charge traps throughout the arra y immediately be- 
fore slewing to the science t arget (e.g. JSeager fc Demingl 
120091: iKnutson et al.ll2009cT ). These data are particularly 
useful for characterizing the ramp effect because they 
allow us to bin pixels according to their median illumi- 
nation level without worrying about added variability in- 
troduced by shifts in the telescope pointing. For a bright, 
compact source such as a star, it is almost impossible to 
disentangle these two effects at the pixel level. 

We divide the data from the preflash images into 
twelve bins, with illumination levels ranging from 1000 — 
4500 MJy Sr -1 . Although these illumination levels are 
generally higher than the peak flux of approximately 
1700 MJy Sr _1 reached in the central pixel of our 
8 /tm HD 189733b observations, we expect that lower 
illumination pixels will display qualitatively similar be- 
haviors on longer time scales. We plot the resulting light 
curves for each bin in Fig. [TT] Although the majority of 
pixels exhibit the typical asymptote associated with this 
effect, we note that the highest illumination pixels dis- 
play an additional behavior which we term "overshoot- 
ing" , in which the flux increases past its equilibrium value 
and then gradually decreases until it reaches its nominal 
steady state. It is possible that this effect is present in 
the lower-illumination light curves as well, but is hidden 
by the larger size of the i nitial ramp. There is s ome ev- 
idence for this theory, as lLaughlin et al.l (|2009[ ) report 
an effect resembling overshooting in their 8 /tm observa- 
tions of th e stellar binary c ompo nent in the HD 80606 
system. In Crossfie ld et all ()2012f ). we also find evidence 
for a similar effect in the 24 /im MIPS array. 

We note that exponential f its to the 8 /tm ramp, such 
as those used in lAgol et all (|2010l) , do not take this ef- 
fect into account. Although we used a more general 
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Figure 11. Measured flux as a function of time in the 8 fim array 
for pixels with illumination levels ranging from 1000 — 2000 MJy 
Sr~ x (bottom curve) to 4000 - 45000 MJy Sr" 1 (top curve), with 
best-fit fourth-order polynomial functions of time (red lines) over- 
plotted for comparison. Although the detector ramp typically re- 
sembles an asymptote that converges on time scales related to the 
median illumination level, the most highly illuminated pixels ex- 
hibit an additional over-shooting effect that could lead to residuals 
in the corrected light curves if not taken into account in our choice 
of functional form to describe the ramp. 

seventh-orde r polynomial f u nction of time to describe 
the ramp in IKnutson et al.l (|2007t ). we only corrected a 
subset of the lower-illumination pixels inside our photo- 
metric aperture. If the uncorrected higher-illumination 
pixels in our aperture displayed this behavior, or if it 
was inadequately described by our polynomial fits to 
the lower-illumination pixels, then we would expect to 
see a local minimum near the start of our observations. 
Based on these data, we conclude that the flux mini- 
mum observed in our corrected data could reasonably 
be attributed to this overshooting effect, rather than the 
planet's phase curve. For the purposes of comparing our 
8 fim phase curve to other w avelengths, we adopt the 
convention of lAgol et al.l (|2010D and trim the first part of 
the light curve where the ramp correction is largest and 
our conclusions correspondingly uncertain. This trim- 
ming does not affect our estimate of the location of the 
flux maximum, but it does prevent us from determining 
the location of the flux minimum in these data. 

4.4. Comparison to General Circulation Models 

In this section we combine our new 3.6 and 
4.5 /im pha se curves with our prev ious observations at 8.0 
and 24 /im (|Knutson et al.ll2009al) and compare these re- 
sults toJhe_p£e^ictMnsof3D general circulation models 
from lShowman et al.l (|2009t) . These models use the MIT- 
gcm to solve the primitive equations in three dimensions, 
and are coupled to a non-grey radiative transfer scheme. 
Showman et al. considered four complementary cases for 
the planet, corresponding to: 1) a solar metallicity atmo- 
sphere and a 1:1 ratio for the planet's rotation and orbital 
periods, 2) a 5x solar metallicity atmosphere and 1:1 ro- 
tation, 3) a solar metallicity atmosphere and a 1:2 ro- 
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Figure 12. Measured phase variations in the 3.6 (blue circles), 4.5 
(green circles), 8.0 (orange open circles from Knutson et al. 2007, 
2009a, orange filled circles from Agol et al. 2010), and 24 fim (red 
diamonds) Spitzer bands after correcting for stellar flux variations 
and instrument effects. Overplotted lines indicate the predictions 
of general circu l ation models for this planet in each band from 
IShowman et~aTI ( 2009), assuming either solar (solid) or 5x solar 
(dashed) metallicity. 
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Figure 13. Same data as in Fig. 1121 The overplotted lines indi- 
cate the predicti ons of general cir culation models for this planet in 
each band from [Showman et al. (2009), assuming either a slowly 
rotating (0.5 X orbital period; solid line) or rapidly rotating (2x 
orbital period; dashed line) planet. 



tation period (i.e., more slowly rotating than the tidally 
locked case), and 4) a solar metallicity atmosphere and a 
2:1 rotation period (i.e., the rapid rotation case). We cal- 
culate phase curves in the 3.6, 4.5, 8.0, and 24 /im Spitzer 
bands for each mode l following the method described in 
iFortnev etHTl (|2006n . We use an interpolated PHOENIX 
model to calculate planet-star flux ratios and confirm 
that this model gives predictions consistent with those 
of the ATLAS model used with the ID models in £14.2.11 
We compare the resulting predictions to our measured 
phase curves in Fig. and Q~3] see Table [H for a more 
quantitative comparison. 

In examining the locations of the local maxima and 
minima, we find that the slow rotator model is the only 
one that provides an acceptable fit (< 3cr disagreement) 
to the measured offsets in the 4.5, 8.0, and 24 /im bands. 
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Table 2 

Phase Curve Comparison 



Source 


3.6 fim 


4.5 /im 


8.0 (im 


24 fim 


A/f /~iTiiy~ii I'm P) nn ft - nv Tt'litT Pfftin^ 

1 VI LLJj VI I bU/1 I b 1 LLLt t>\^L vU/l 1 L LLJb 1 LLLLt/U 










Measured 13 


0.1466% ±0.0040% 


0.1787% ± 0.0038% 


0.344% ± 0.004% 


0.598% ± 0.038% 


lx Solar Model 


0.1397% 


0.1948% 


0.328% 


0.552% 


5x Solar Model 


0.1563% 


0.2170% 


0.363% 


0.596% 


Slow Rot. Model 


0.1326% 


0.1911% 


0.326% 


0.553% 


Fast Rot. Model 


0.1270% 


0.1806% 


0.310% 


0.532% 


TlniiQin/^ JPliiT P irti n /PnnQp — ^ )^ 

J_y LL U 'J t/LLC 1 b UjiJb H.LLU b\J I i / bUrJ O \J *\J I 










lx Solar Model 


0.1194% 


0.1760% 


0.301% 


0.523% 


5x Solar Model 


0.1423% 


0.2035% 


0.347% 


0.581% 


Slow Rot. Model 


0.1305% 


0.1884% 


0.321% 


0.547% 


Fast Rot. Model 


0.1092% 


0.1656% 


0.285% 


0.507% 


A/TtTit fTi 1 1 TT1 P] f~l pi flT* PI 1 ) 1 *Y P n ~fl 
1V1 LI III 1 bill 1 L 1 LLLlbCtj tJLLLl i t LLJ, 1 LtLobU 










Measured 15 


0.0224% ± 0.0073% 


0.0805% ± 0.0097% 


< 0.311% c 


0.443% ± 0.027% 


lx Solar Model 


0.0384% 


0.1063% 


0.158% 


0.375% 


5x Solar Model 


0.0382% 


0.0912% 


0.150% 


0.359% 


Slow Rot. Model 


0.0413% 


0.1124% 


0.173% 


0.399% 


Fast Rot. Model 


0.0535% 


0.1248% 


0.189% 


0.414% 


Nigntside tlux Katio (Fnase= 0/ 










lx Solar Model 


0.0605% 


0.1263% 


0.194% 


0.413% 


5x Solar Model 


0.0557% 


0.1073% 


0.175% 


0.385% 


Slow Rot. Model 


0.0424% 


0.1128% 


0.174% 


0.400% 


Fast Rot. Model 


0.0800% 


0.1423% 


0.229% 


0.450% 


Minimum Flux Offset [hj 










Measured 


-6.43 ±0.82 


-1.37 ± 1.00 


< 10 


< 


lx Solar Model 


-8.4 


-8.9 


-7.8 


-7.9 


5x Solar Model 


-6.9 


-6.9 


-5.9 


-5.9 


Slow Rot. Model 


-2.7 


-2.3 


-2.4 


-2.3 


Fast Rot. Model 


-10.6 


-10.6 


-10.1 


-10.1 


Maximum Flux Offset [h] 










Measured 


-5.29 ±0.59 


-2.98 ±0.82 


-3.5 ±0.4 


-5.5 ± 1.2 


lx Solar Model 


-7.8 


-7.6 


-7.3 


-7.2 


5x Solar Model 


-6.0 


-5.2 


-4.4 


-4.0 


Slow Rot. Model 


-2.5 


-2.7 


-2.8 


-3.1 


Fast Rot. Model 


-9.0 


-8.6 


-7.9 


-7.6 



a Because our maximum fluxes are located close to a phase of 0.5, we approximate the maximum fluxes as the measured secondary eclipse 
depths. The minimum flux is then simply the maximum flux minus the measured phase curve amplitude. Model fluxes are calculated at a 
single time corresponding to cither the average location of the flux maximum/minimum or an orbital phase of cither (night side) or 0.5 
(day side). 

Fluxes are estimated after removing fitted (3.6 and 4.5 /im) or predicted (8.0 and 24 /im) stellar flux variations due to spots. 
c la upper limit based on the minimum flux estimate of 0.261% ± 0.025% from lKnutson et al.l J2009al1 . 



The 3.6 (im band displays a larger offset that is incon- 
sistent with this model at the 5a level, and is instead 
best-matched by the offset in the 5x solar model. We 
note that there are a number of model parameters in ad- 
dition to rotation rate that could reduce the offsets of 
the maxima and minima given the complex interplay of 
radiative, chemical, and advective timescales. For ex- 
ample, increasing atmospheric drag in the region near 
the photosphere would decrease the wind speeds and re- 
duce t he corresponding offsets (e.g.. iRauscher fc Menoul 
l2012al) . By controlling the opacities and therefore the 
photosphere pressures that contribute flux in a given 
bandpass, atmospheric chemistry also plays a large role 
in determining the offsets of the observed maxima and 
minima. Our models assume equilibrium chemistry, but 
it is likely that disequilibrium pro cesses are important 
for HD 189733b (IMoses et al.H2011l) . 

Although the slow rotator model provides a reasonably 
good match to the measured locations of the flux max- 
ima and minima, we must next consider whether or not 



it accurately predicts the planet's emission spectrum at 
these two points in time. We find that the our measured 
maximum fluxes differ by [±11%, -7%, ±5%, ±8%] from 
their predicted values in the [3.6, 4.5, 8.0, 24] bands, 
with a significance of [3.5er, — 3.3er, 4.5c, 1.2c], respec- 
tively. On the planet's night side, the minimum fluxes 
differ by [-46%, -28%, ±11%] from their predicted val- 
ues in the [3.6, 4.5, 24] /im bands, with a corresponding 
significance of [— 2.7<r, — 3.3<r, 1.6c], respectively. The 
2tT upper limit on the 8 /im minimum is consistent with 
all models. We show the predicted spectra at times of 
minimum and maximum flux in Fig. 1 141 and list the cor- 
responding planet-star flux ratios in the Spitzer bands in 
Tabled 

We next consider whether or not non-equilibrium 
chemistry might be able to explain the disagreement be- 
tween the measured and predicted planet-star flux ratios. 
The planet's 3.6 and 4.5 /im nightside fluxes display the 
single largest discrepancies, with valuees approximately 
one third to one half those predicted by our preferred 
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Figure 14. Comparison of measured broadband planet-star flux 
ratios (filled bla ck circles) at 3.6 and 4.5 jim (this paper), 5.6 , 
16, and 24 nm llCharbonneau et aJ1l20 08; Knut son etaT] [2009a). 
and 8 /im I IKnutson et al.H2009al ; lAeol et al.112010) , aiid an IRS 
spectrum (open grey circles, Grillmair et al. 2008) vs. 3D model 
emission spectra at times of minimum and maximum flux, approx- 
imately corresponding to the night and day sides of the planet. 
Solid lines indicate the model predictions for the lx solar (blue), 
5x solar (green), slow rotator (orange) and fast rotator (red) cases, 
with the bandpass-integrated fluxes for each model overplotted as 
filled squares with the same color. 

model. As discussed in lCooper fc Showman! (|2006| ). the 
time scale for converting between CO and methane is 
long relative to the vertical mixing time scales in these 
atmospheres, and as a result the relative abundances of 
these two molecules should reflect their values in the 
planet's deeper atmosphere (~ 1 — 10 bars) where chemi- 
cal quenching can occur. For planets in the temperature 
regime of HD 189733b, Cooper & Showman's models sug- 
gest that approximately 20% of the total carbon should 
exist in the form of methane on both the day and night 
sides, with the remainder residing primarily in CO. If we 
compare this to chemical equilibrium predictions for HD 
189733b's night side, we find that this leads to an ex- 
cess of CO and a deficit of CH4 at the pressures probed 
by our observations. Increasing the amount of CO in our 
nightside model would suppress the 4.5 /Ltm flux, bringing 
it into good agreement with our measured value in this 
band. This result is also consistent with our constraints 
on the planet's transmission spectrum (see §4.2j) . which 
suggest the presence of additional CO absorption at the 
day-night terminator. Although the 3.6 /jm nightside 
flux is also lower than predicted, we note that the slow 
rotator model provides a poor match to the measured lo- 
cations of the flux maximum and minimum in this band, 
and it is therefore not surprising that it would also fail 
to match the measured flux on the planet's night side. 

Although non-equilibrium chemistry offers a possible 
solution for HD 189733b's ni ght side, the pictu re on 
the day side is more m ixed. iMoses et al.l ([201 1[ ) and 
iVisscher fc Mosesl ([201 ll ) find that their non-equilibrium 



models for this planet's day side have more methane and 
less CO than the equilibrium predictions, as expected in 
the limit of strong vertical mixing (K zz — 10 11 ). How- 
ever, the net effect of this non-equilibrium chemistry in 
their models was to decrease the dayside fluxes in the 3.6, 
8.0, and 24 /jm bands, while preserving a nearly identi- 
cal flux in the 4.5 /nm band; unfortunately this is the 
opposite of the change needed in this case. 

5. CONCLUSIONS 

In this paper we present new full-orbit and near- 
continuous phase curve observations of the hot Jupiter 
HD 189733b in the 3.6 and 4.5 fim Spitzer bands, which 
allow us to characterize the atmospheric circulation pat- 
terns and corresponding chemistry of HD 189733b's at- 
mosphere. These data include one transit and two sec- 
ondary eclipses in each bandpass, which we use to derive 
an improved estimate for the planet's orbital ephemeris 
and wavelength-dependent transmission and emission 
spectrum in these bands. We confirm that the planet's 
4.5 /im transit depth is 3<r smaller than at 3.6 /itm, 
consistent with the presence of excess CO at the day- 
night terminator, a lthough our precisio n is comparable 
to that reported bv lDesert et al.1 fj2009f ) . We obtain im- 
proved estimates for HD 189733b's dayside planet-star 
flux ratio of 0.1466% ± 0.0040% in the 3.6 /im band and 
0.1787%±0.0038% in the 4.5 /zm band, corresponding to 
brightness temperatures of 1328±11 K and 1192±9 K, re- 
spectively; these are the most accurate secondary eclipse 
depths obtained to date for an extrasolar planet. Our 
new 3.6 /nm secondary eclipse depth is 7.5g s malle r than 
the value reported in iCharbonneau et ail ([2008), but 
we find that the uncertainties in this previous measure- 
ment, which assumed Gaussian noise, are likely underes- 
timated. We conclude that there is no convincing ev- 
idence for time variability in the measured secondary 
eclipse depths or times, consist ent with the upper lim- 
its derived by lAgol et al.l (|2010l ) using a more extensive 
8 /zm data set. 

We combine our new 3.6 and 4.5 /im phase curves with 
previously published observations at 8.0 and 24 fim in 
order to characterize the planet's emission spectrum as 
a function of orbital phase. We find that the times of 
minimum and maximum flux occur several hours earlier 
than predicted for an atmosphere in radiative equilib- 
rium, consistent with the eastward advection of gas by 
an equatorial super-rotating jet. The locations of the 
flux minima in our new data differ from our previous 
observations at 8 /im, and we present new evidence in- 
dicating that the flux minimum observed in the 8 /im is 
likely caused by an over-shooting effect in the 8 /im array. 
We fit the planet's day side spectrum with a ID radia- 
tive transfer model from lBurrows et all ([20081 ) where the 
amount of energy transported to the night side is left as a 
free parameter and find that the corresponding nightside 
spectrum is in good agreement with our measured val- 
ues. This serves to validate studies (e.g.. iCowan fc Aeoll 
120 111 ) that seek to constrain circulation patterns on hot 
Jupiters based on dayside spectra alone. 

We then compare our phase curves to the predictions 
of 3D general circulation models from [Showma n et al.l 
(2009) and find that we require models with either a 
slower-than-synchronous rotation rate or increased drag 
at the bottom of the atmosphere in order to match the 
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small measured offsets in the locations of the phase curve 
maxima and minima at 4.5, 8.0, and 24 /im. We also 
find that HD 189733b's 4.5 ^m nightside flux is 3.2a 
smaller than predicted by these models, which assume 
that the chemistry is in local thermal equilibrium. We 
conclude that this discrepancy is best-explained by ver- 
tical mixing, which should lead to an excess of CO and 
correspondingly enhanced 4.5 /im absorption in this re- 
gion. This result is consistent with our constraints on 
the planet's transmission spectrum, which also suggest 
excess absorption in the 4.5 ^m band at the day-night 
terminator. 

Looking ahead, it is clear that the questions re- 
garding atmospheric chemistry and circulation patterns 
on tidally locked planets will continue to recur as 
the field shifts towards studies of smaller and more 
earth-like worlds. Current studies of the atmospheres 
of low-mass planets such a s GJ 1214b (iBean et al.1 
l20~LOl I201H: IBerta et all 12011 IDesert et al.ll2011bT) and 



GJ 436b fe.g.. iStevenson et al. I 120101: iBeaulieu et al.l 



I201H Knutson et all 1201 lh focus almost exclusively on 
systems with M star primaries, as the lower stellar effec- 
tive temperature and smaller stellar radius result in pro- 
portionally larger transit and secondary eclipse depths. 
For late M stars, the location of the habitable zone is 
within the re gion in which we wo uld expect tidal lock- 
ing to occur (|Kasting et al.l fl993h ; it is therefore likely 
that the first atmosphere studies of potentially habitable 
worlds with the James Webb Space Telescope will focus 
on these tidally locked systems. 
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